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ABSTRACT 



Context. To derive the history of star formation in the Universe a set of calibrated star formation rate tracers at different wavelengths 
is required. The calibration has to consistently take into account the effects of extinction, star formation regime (short or long-lived) 
and the evolutionary state to avoid biases at different redshift ranges. 

Aims. We use evolutionary synthesis models optimized for intense episodes of star formation to compute a consistent calibration of 
the most usual star formation rate tracers at different energy ranges, from X-ray to radio luminosities. 

Methods. We have computed the predicted evolution of the different estimators taking into account nearly-instantaneous and continu- 
ous star formation regimes and the effect of interstellar extinction (attenuation at high energies, thermal reradiation in the far infrared). 
We have also considered the effect of metallicity on the calibration of the different estimators. 

Results. A consistent calibration of a complete set of star formation rate tracers is presented, computed for the most usual star-forming 
regions conditions in terms of evolutionary state, star formation regime, interstellar extinction and initial mass function. We discuss 
the validity of the different tracers in different star formation scenarios and compare our predictions with previous calibrations of 
general use. 

Conclusions. In order to measure the intensity of star formation episodes we should distinguish between nearly-instantaneous and 
continuous star formation regimes. While the star formation strength (M Q ) should be used for the former, the more common star 
formation rate (M Q yr -1 ) is only valid for episodes forming stars at a constant rate during extended periods of time. Moreover, even 
for extended star formation episodes, the evolutionary state should be taken into account, since most SFR tracers stabilize only after 
100 Myr of evolution. 



Key words. Stars: formation - Galaxies: evolution - Galaxies: fundamental parameters 
function - Galaxies: starburst - Galaxies: stellar content 
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1. Introduction 

The availability of a continuously increasing observational data 
set on galaxies at higher and higher redshifts allows the study 
of the history of star formation at cosmic scales, approaching 
already (presently up to z > 7) the ages of reionization when 
the first population of massive stars started to ignite. The star 
formation rate (SFR) of these galaxies at different ages is de- 
rived from the so-called star formation rate tracers, spatially in- 
tegrated parameters at different wavelengths that allow the deter- 
mine of the strength of the star formation episode by c omparison 
with th e predictions of evolutionary synthesis models. Ke nnicutil 
(1998) discussed the uncertainties inherent to the use of opti- 
cal broad-band colors, and presented an updated calibration of 
the UV continuum, recombination lines (or number of ionizing 
photons emitted per unit time), forbidden lines and total far in- 
frared emission that has become de facto the reference standard 
calibration of these estimators. 

While widely used, i t is often n e glecte d that the SFR tracers 
calibration proposed by Ke nnicutil ( 1998b was computed for a 
very specific star formation history, i.e., a star formation episode 
running at a nearly constant rate during a long enough period of 
time (> 100 Myr), so that an equilibrium is reached between the 
number of massive stars dying and igniting. Under these con- 
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ditions most of the parameters considered reach an equilibrium 
state with almost no evolution with time. This continuous burst 
scenario is certainly valid for many galaxies, especially large spi- 
rals in which massive star formation takes place at spatially inte- 
grated nearly constant rates during long periods of time (though 
at different locations within their spiral arms). But star forma- 
tion seems to proceed in a more bursty regime in other galaxies, 
showing extremely high present-day star formation rates, indi- 
cating that they might be experiencing a nearly instantaneous 
(i.e., extended over few million years) massive star formation 
episode. Under these circumstances, the derivation of star for- 
mation rates by comparison with the predictions by continu- 
ous bursts might be erroneous. Moreover, even the concept of 
a star formation rate, measured as M yr _1 , might be mislead- 
ing in these cases. The intensity of these massive star formation 
episodes should be better parameterized by the star formation 
strength (SFS), measured directly in units of M Q , indicating the 
total amount of mass having been transformed into stars. 

On the other hand, as noted by [Rosa-Go nzalez et "ai1 d2002l) . 
different studies have indicated the significant role played by 
dust in the estimates of star formation rates. Extinction by dust 
significantly alters the integrated multiwavelength spectrum es- 
caping from a star-forming region, weakening the far-UV con- 
tinuum, but also boosting the emission in the far infrared range. 
It is critical therefore to take into account the effects of extinction 
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by dust in a consistent way, since the calibration of a multiwave- 
length set of SFR and SFS tracers will be a strong function of the 
dust abundance. 

In this paper we present an updated calibration of a com- 
plete set of star formation rate tracers at different wavelengths, 
from X-rays to radio, computed in a consistent way with state- 
of-the-art evolutionary synthesis models considering different 
star formation histories, and taking into account the effects of 
dust extinction at all wavelengths. In Sect. 2 we summarize the 
properties of the tracers considered in this work, and in Sect. 3 
we describe the evolutionary synthesis models we have used. 
In Sect. 4 we present the calibration of the different estimators, 
discuss their sensitivity to the properties of the star formation 
scenario and compare our predictions with other evolutionary 
codes. Finally, Sect. 5 summarizes our main conclusions. 

2. Star formation rate tracers 

We have selected a large sample of SFR and SFS tracers cover- 
ing a broad wavelength range, from X-rays to radio, and associ- 
ated to different physical processes directly linked to the strength 
of the star formation episode. We want to remark that most of 
these tracers may be contaminated by emission not related to the 
present episode of star formation, but associated to an existing 
AGN, underlying older stellar populations, evolved but hot, low- 
mass stars, etc. During this work we will consider that the dif- 
ferent star formation rate tracers have been corrected by the user 
before applying the calibrations. Since these corrections will not 
always be accurate or even possible, the SFR values so derived 
have to be taken with the corresponding caution. The reliabil- 
ity of the estimates will be larger the more SFR tracers can be 
used simultaneously for a given object, since the contaminating 
effects have different strengths at different energy ranges. 

2. 1 . Far infrared luminosity (Lms.) 

Most radiation emitted by the young, massive stars is absorbed 
by the dust particles surrounding the star formation region, prin- 
cipally in the UV. The heated dust reemits this energy in the far 
infrared (FIR) range (which we will consider as 1 - 1000 fim), 
especially within the 10 - 120 //m domain. The fraction of en- 
ergy emitted by these stars in the UV/optical which is absorbed 
by dust and reemitt ed as FIR radiation is close to 1 for values of 
E(B-V) above 0.5 dMas-Hesse & Kunthlfl99ll) . Hence, the FIR 
luminosity is in general a good estimator of the bolometric lu- 
minosity of the massive stars, and thus SFR(FIR ) is on e of the 
most used and reliable SFR estimators. Ken nicuttld 19981) derived 
a calibration of Lffr using population synthesis models of ex- 
tended bursts in their asymptotic phase (when the rate of mas- 
sive star formation and death is balanced), assuming the opti- 
cally thick case, but stating that the expression strictly applies to 
bursts younger than 100 Myr. This calibration would be altered 
if corrections for the dust heating from old stars in quiescent 
galaxies were made. 

Besides the total FIR luminosity, some other calibrations 
make use of the SED (spe ctral energy distribution) of t he burst 
at certain FIR wavelengths. lAlonso-Herrero et al.l d2006l) derived 
an empirical, non-linear calibration for S FR(L(24 fim)) by com- 
bining L(24 //m) observations of Luminous InfraRed Galaxies 
(LIRGs), Ultraluminous InfraRed Galaxies (ULIRGs), normal 
galaxies and HII regions within M51, and t he SFR(Pa) val- 
ues derived from the SFR(Ha) cali bration from | Kennicuttl (fl 998 ) 
(assuming case B recombination). ICalzetti et all (120071) found a 
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Fig.l. Evolution of LW(WRbump) predicted by CMHK02 
for IB (solid lines) and EB (dashed lines) models both for 
Z = 0.020 (top) and Z = 0.008 (bottom), together with 
the histogram of the LWfWRbump ) data of the sample from 
iBrinch mann. Kunth&Durrel (120081) . 



similar relation based on luminosity surface density values of 
HII knots from 33 nearby galaxies. 

Lfir can become contaminated by the presence of an AGN, 
especially when dealing with spatially-integrated luminosities of 
distant, unresolved galaxies. Depending on the redshift of the 
galaxy, the width of the Balmer or Lyman emission lines should 
be checked to exclude AGN dominated objects. In galaxies dom- 
inated by massive star-forming episodes, the SFR values derived 
from Lfir and from Ljjy, computed with the same E(B-V) levels, 
should be consistent. 

2.2. Ionizing power (N Lyc ) and emission lines 

Massive stars are conspicuous sources of ionizing photons, and 
therefore their number A^yc could trace the burst intensity di- 
rectly. However, practically all ionizing photons are absorbed by 
gas and dust within the nebular region. Recombination of the 
hydrogen atoms ionized by this radiation and free electrons pro- 
duces intense emission lines, such as Ha or HB, commonly used 
as SFR estimators since th eir luminosity is p r oportional t o A^lvc 
under basic assumptions (Osterbrock 1989|. IKennicutl {l998) 
pr ovided an expression fo r SF R(Nlvc) based on the calibrations 
by iKennicutt et al.l d 19941) and iMadau et al.l dl998l) . and used it 
to derive a calibration of SFR(Ha) assuming case B recombi- 
nation and T P = 10 4 K. A som e what different expression was 
derived by Rosa-Gonzalez et al. (2002). They studied a sample 
of HII, starburst and blue compact galaxies, in order to obtain 
empi rical SFR calibrations based on L(Ha), L([OII]3727) and 
Luv Rosa-Gonzal ez et all d2002i) obta ined an SFR(Ha) coef- 
ficient 40% lower than the value from iKennicuttl d 19981) since 
it was based on the observed Ha luminosity, not corrected for 
Balmer absorption due to the stellar population. 

Ha has the disadvantage of falling outside the opti- 
cal spectrum when studying medium-redshift sources (z > 
0.4). Studying a samp le of 412 local star-forming galaxies, 
Moustakas et al ] (I2004 showed that, although weaker and more 
affected by extinction, HB line can a lso be used to esti mate SFR. 
Applying the L(Ha) calibration by IKennicutl {1998), they cal- 
culated the SFR values using both L(Ha) and L(Hy8), assuming 
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case B recombination (i.e. L(Ha)/L(HjS) = 2.86). They obtained 
very similar results, with a 40% scatter. 

The luminosity of the Lya line can also be used as a star 
formation rate estimator, since it is in principle a linear function 
of A^Lyc, with the advantage that it is visible in the optical range 
for galaxies at a redshift of z > 2.0. Nevertheless, Lya photons 
are affected by resonant scattering in neutral Hydrogen, strongly 
affecting the intensity of the line, as we will discuss later. 

2.3. Continuum luminosity 

Most radiation from young, massive stars is emitted in the UV 
range (912 to around 3000 A) and could be used as a reliable SFR 
and SFS estimator since the contamination from older stellar 
populations is low or even negligible below 1800 A. However, 
attenuation correction plays an importan t role here, since UV 
photons are severly affected by extinctio n. lKennicuttl (fl998) uses 
the calibration from Madau et al. (1998) in order to obtain an ex- 
pression for SFR as a function of L Y , assuming an extended star 
formation process with a n age larg e r than 100 Myr and solar 
metallicity. According to iKennicuttl dl998l) . the expression can 
be used in the wide range of 1500 - 2500 A, where the spec- 
trum of a stellar population with a Salpeter initial mass function 
(IMF) is rather flat (as a function of frequency). 

Fitting GALEX and SDSS (Sloan Digital Sky Survey) data of 
nearly 50000 galaxies w ith SEDs construct e d with the popula- 
tion sy nthesis models of [Bruzual & Charlol d2003l) . ISalim et all 
(2007) derived an expression for SFR(\JV) which yields SFR 
values lower by 30% than those obtained with the IKennicuttl 
(1998) calibration. They claim that the reasons for this differ- 
ence are the low metallicity of their sample (Z = 0.016), the 
different star formation histories of th e objects of the sample and 
some intrinsic d i fferen ces between Bruz ual & Charloil {2003) 
and lMadau et all dl998l) models 

Some authors ( Buat et al. 19991 : iHirashita. Buat & Inoud 
l2003t l!d esias-Paramo et aUl006) have opted for merging in 
a single expression the observed values of UV and FIR lumi- 
nosities. The latter component would represent the main bulge 
emission of massive stars, while the former would correct for 
the radiation which eventually escapes from the burst before it is 
absorbed by dust. This way, prior assumptions about extinction 
would not be necessary, but both the FIR and the UV luminosi- 
ties would be needed, so that the method is not easily applicable 
for high-redshift galaxies. 

IKennicuttl 0998) concluded that broad-band luminosities in 
the optical are poor SFR estimators, since the optical continuum 
is contributed by stars at very different evolutionary states and 
affected by potentially different degrees of extinction. This is 
especially true for large spiral galaxies with a complex history 
of star formation. Even the UV continuum might be contami- 
nated by blue horizontal-branch stars in those objects dominated 
by an old underlying stellar population, specially in the case of 
early-type galaxies. Nevertheless, continuum luminosities in the 
optical-near infrared bands could be reliable estimators for com- 
pact objects whose continuum is dominated by the present burst 
of star formation. We have calibrated SFR and SFS with the lu- 
minosities at 1500 A, 2000 A, 3500 A (U), 4400 A (B), 5500 A 
(V) and 2.2 /urn (K). 

2.4. X-ray luminosity (L x ) 

Mechanical energy released by the massive stars' stellar winds 
and supernova explosions heats the surrounding gas, origi- 



nating a diffuse X-ray emission which peaks in the soft X- 
ray range. This diffuse emission adds to the X-ray point 
source radiation from X-ray binaries and supernova remnants 
(SNR). Whereas the emission from X-ray binaries dominates 
the hard X-rays range (2-10 keV), emission from dif- 
fuse hot gas drives the radiation in the soft range (0.2 - 2 
keV) dCervifio. Mas-Hesse & Kunthll2002l) . Many authors have 
claimed recently that the X-ray emission from star-forming re- 
gions should be a direct function of the burst intensity, and sev- 
eral SFR estimators based on different components and energy 
ranges of this radiation have been published. Analyzing the total 
X-ray emission from local starburst galaxies extracted from the 
iHo et al.l d!997[) atlas, which had been observed by ASCA and/or 
BeppoSAX, iRanalli. Comastri & SettH d2003l) found both soft 
and hard X-rays SFR linear express ions using FIR a s a proxy, 
based on the SFR(FIR) re lation b y Kennicutt (1998). A simi- 
lar study was performed bv lTullmann et al.l (2006) using XMM- 
Newton and Chandra data, but in t his case obvious poi nt sources 
were removed. On the other hand. iGrimm et al] d2003l) used the 
hard emission from X-ray binaries of Chandra-resolved, nearby 
late-type/starburst galaxies, t o derive a relation wh i ch bec omes 
non-linear at low SFR values. iRosa-Gonzalez et all d2009) have 
recently analyzed XMM-Newton observations of a sample of 14 
star-forming galaxies with high SFR (0.2 - 160 M Q yr _1 ), as de- 
rived from radio, FIR and UV tracers, confirming that the SFR 
derived from soft X-rays is comparable to that obtained from Ha 
luminosities. 

Using evolutionary synth esis models, 

iMas-Hesse. Oti-Floranes & Cervind d2008l) proved the im- 
portance of the evolutionary state of the burst when trying to 
ascertain the SFR value via the soft X-ray luminosity, since 
it can increase between half and one order of magnitude 
within 20 Myr in the models with constant SFR. They de- 
rived two expressions, both for young (~ 10 Myr) and more 
evolved extended star formation bursts (~ 30 Myr). They 
also provided a star formation strength expression for nearly 
instantaneous star-forming episodes, yielding the initial, total 
stellar mass of gas transformed into stars. We have included 
the calibration of Lx in this paper for c omplet eness. We refer 
to IMas-Hesse. Otf-Floranes & Cervifiol |2008) for a deeper 
discussion. 



2.5. Radio luminosity (L mA ) 

A very tight correlation between Lfir and radio emission has 
been found over the past decades, common to a great variety 
of objects, such as starburst galaxies, normal spiral galaxies, 
blue compact dwarfs ( BCDs), E/SO gala x ies, irregular galax- 
ies, H II regions, etc... dMas-Hessd [1991 ICondonT[T992T iBeiH 
2003), which proves that radio luminosity might be a reliable 
star formation rate estimator. Radio luminosity in star-forming 
regions is composed of both a thermal (L tn ) and a non-thermal 
(L ntn ) component. The former appears as a byproduct of free- 
free interactions, i.e. Bremmstrahlung, and free-bound transi- 
tions between the constituents of the ionized g as and its value 
can be expressed in terms of A^ yc (Rubin 1968), while the non- 
thermal emission is synchrotron in nature, emitted by electrons 
after being accelerated by supernovae explosions, and it is usu- 
ally assumed to be proportional to the rate at which SN explode 
(IUlvestadlll982l) . Both contributions hav e different sp ectral in- 
dex values, which are a tn = -0.1 for L tn ([Rubin 1968) and a ntn 
= -(0.8 - 0.9) for L nt h (see Mas -Hessel (1 1992b and references 
therein). SFR cal i bratio ns based on both L ntn and L tn can be 
found in ICondonl (1 19921) . assuming a - -0.8 for the former 
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dCondon & Yin|[T990l) and an electronic temperature T e - 10 4 
K and no dust abso rption for the la tter. Applying the calibration 
of SFR(FIR) from iKennicuttl (1 1 9981) to a sample with objects of 
very diverse natures (spiral galaxies, starburst galaxies, BCDs, 
irregular galaxies, etc), and after consideri ng the tigh t correla- 
tion between infrared and radio emissions. iBelll (2003) obtained 
a non-linear expression for SFR(L m &) using Lfir as a proxy. 



3. Evolutionary synthesis models 

In order to calibrate the different SFR tracers in a consistent 
way, we have computed their expected values under different 
scenarios using the evolutionar y population synthesis models of 
ICeryifio, Mas-Hesse & Kunthl (2002b (hereafter CMHK02 mod- 
el Jl), [Leitherer et alj d 1999b (S tarburst99, hereafter SB99 mod- 
elo and ISchaererl (12002. 2003) (SC02 models^). 

CMHK02 m odels, which are based on the models 
by lArnault et al.l 



- hVlas-Hesse & Kunthl d 1991 1) and 

ICervino & Mas-Hesse (11994 ), compute the evolution of a young 
population of massive stars which are formed at the same time 
(instantaneous bursts, IB, also referred to in the literature as sin- 
gle stellar populations - SSP) or which form at a constant rate 
during an extended period of time (extended bursts, EB). The 
different observables are calculated for the first 30 Myr after the 
onset of the burst. They are scaled to the mass of gas transformed 
into stars at the start of the burst in IB models (measured in M ), 
and to the mass of gas transfomed into stars per unit time for EB 
models (M yr -1 ). The initial mass function (IMF) of the stellar 
population is defined by a power-law within the mass range of 
2-120 M with the slopes a = 1, 2.35, 3. The models include 
different metallicities, Z/Z© = 0.05, 0.2, 0.4, 1, 2. For this work 
we have considered a nominal scenario with <p(m) ~ m~ 235 for 
the IMF (i.e. a Salpeter IMF), and solar metallicity Z = Z©. The 
presence of binary systems has not been considered. Chemical 
evolution is not treated self-consistently, since all stellar genera- 
tions in EB models have the same metallicity. This should not be 
a problem, since th e effect becomes important at larg e ages (~ 1 
Gyr), as shown bv lFioc. & Rocca-V olmerangel dl997l) . 

SB99 synthesis models (Leith erer et al.ll 1 999b ar e based on 
the original models by ILeitherer & Heckm an (1995), but have 
been continuously updated with new evolutionary tracks and 
stellar atmospheres. Several evolutionary tracks with differ- 
ent metallicity values can be used, and for consistency with 
CMHK02 models, we opted for the Geneva standard mass-loss 
tracks. SB99 allows us to compute the observables up to 1 Gyr 
after the onset of the burst as well as the use of a step-IMF and 
different mass limits. As in CMHK02 models, both EB and IB 
star formation regimes are considered, and metallicity-evolution 
is not taken into account. We have used SB99 to compute the 
predicted values of the different parameters for EB episodes that 
have formed stars during more than 200 Myr at a constant rate, 
since the models by CMHK02 are available only for the first 
30 Myr. We have checked that where they overlap in time, the 
predictions by both sets of models are perfectly consistent. For 
the IB star formation regime we have preferred to use CMHK02 
models, since they are based on Monte Carlo simulations of the 
IMF and consider therefore in a more realistic way the stochas- 
tical nature of the massive stellar populations. Moreover, since 
they provide directly the expected X-ray luminosity, as discussed 



1 Downloadable from http://www.laeff.inta.es/users/mcs/SED/ 

2 Downloadable from http://www.stsci.edu/science/starburst99/ 

3 Downloadable from http://obswww.unige.ch/sfr/sfr_tls/pop32/ 



below, the other SFR tracers should be computed with the same 
models for consistency. 

Finally, SC02 models have been used for completeness to 
compute Ljjv and /V Lyc f° r population III (Pop. Ill) stellar pop- 
ulations, which could correspond to the scenario valid for star- 
forming galaxies at very high redshift, where metallicity is ex- 
tremely low (Z ~ 0). 

While most SFR tracers calibrations are based on the predic- 
tions by evolutionary models which are computed for extended 
star formation processes where gas is transformed into stars at 
a nearly constant rate during long periods of time (hundreds to 
thousands Myr), we want to stress that star formation seems to 
proceed in the form of short, nearly instantaneous bursts in many 
galaxies. It is generally difficult to distinguish a short, but young 
starburst from an evolved, extended star formation process, since 
most observable parameters are related to just the most mas- 
sive and therefore young stars in the region. Nevertheless, a 
long-lasting star formation process produces an accumulation of 
medium-low mass stars with long lifetimes, which can dominate 
the optical continuum. As an effect, these kind of episodes show 
systematically lower equivalent width values of HyS and Wolf- 
Rayet star bumps at around 4686 A (hereafter EW(WRbump )) 
than young instantaneous bursts (Cervino & Mas-Hes se 1994). 

Figure Q] shows that continuous star formation mod- 
els in the asymptotic phase cannot reproduce the high 
wing of the distribution of £W(WRbump ) values mea- 
sured by iBrinchmann. Kunth & Durre _ d2008l) on their sam- 
ple of 570 star-forming galaxies. Instantaneous bursts, on 
the other hand, of course predict t he hi gh ZiW(WRbump) 
range (see also [Mas-Hesse & Kunthl dl999b ). The sample by 
IBrinchmann. Kunth & Durretl d2008b is clearly biased towards 
star-forming galaxies with large Wolf-Rayet bump equiva- 
lent widths, i.e., experiencing short star formation episodes, 
but other authors have found evidences of short-lived bursts. 
iPellerin & Robert! (2007), for example, find that EB models can 
not reproduce the FUSE far-UV data of the 24 starburst galaxies 
of their sample, whereas there is a much better agreement with 
the predictions by IB models. 

We conclude therefore that star formation proceeds in the 
form of nearly instantaneous, short bursts in a significant frac- 
tion of objects. In the case of very large, non-resolved galaxies, 
the superposition of individually instantaneous bursts at different 
evolutionary states can mimic the properties of an extended star 
formation process at a nearly constant rate, while large objects at 
high redshift could sustain intense star formation processes dur- 
ing truly extended periods of time. When deriving the intensity 
of the different star formation episodes in the history of a galaxy, 
its star formation regime has to be carefully taken into account. 
In the case of a nearly instantaneous burst, the ionizing power 
decreases very fast with time, and becom es negligible after the 
first 8-10 Myr (Ma s-Hesse & Kunth|[l99lb . Blindly applying an 
SFR estimator related to Ni_ yc , calibrated for an extended star 
formation regime, can lead to discrepancies of several orders of 
magnitude depending on the evolutionary state. Moreover, even 
the concept of a star formation rate as such becomes completely 
misleading in these cases, since formation of massive stars might 
have ended already a few million years ago. The concept of a 
star formation strength, SFS, measured as the total mass of gas 
transformed into stars during the burst, is much better adapted to 
this scenario. Unfortunately, deriving the SFS by comparing the 
observations with the predictions of synthesis models requires a 
preliminary estimate of the evolutionary state of the burst, and 
this can be obtained only when parameters like the equivalent 
width of H/3 can be measured. We have therefore calibrated our 
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set of SFR and SFS tracers for both star formation regimes, ex- 
tended and instantaneous respectively, and for different evolu- 
tionary states of the burst. 

Formation of massive stars is an inherently stochastic pro- 
cess, so that when a relatively low number of stars are formed, 
we can not g uaran tee that the whole IMF will be filled. 
ICervino et alj d2002l) showed that stochastic effects become 
small enough only for starbursts that have transfomed more than 
10 5 " 6 M Q of gas into stars. The star formation estimators should 
therefore not be applied to smaller star-forming regions, where 
the stochastic effects could significantly affect the integrated 
emission. 

In the sections below we describe the way the different es- 
timators have been compu t ed. More details can be fou nd in 
Mas-Hesse & Kunthl d 19911). ICervino & Mas-Hessel dl994l) and 
Cervifio. Mas-Hesse & Kunthl d2002T) ~ 



3. 1 . Far infrared luminosity (Lm.) 

The FIR emission predicted by CMHK02 models is calculated 
with the assumption that interstellar dust remains in thermal 
equilibrium. All energy absorbed, mainly stellar UV contin- 
uum, is therefore re radiated in the FIR range (1 - 1000 fj.no). 
ICalzetti et al.l (2000) noticed that this parametrization is higher 
by up to 75% than the value obtained using the expression by 
iHelou et ail dl988l) . which models the IR output within 42.5 - 
122.5 fim via the IRAS fluxes at 60 fim and 100 fim, assum- 
ing a single dust c omponent with T - 2 - 80 K. The Galactic 
extinction law by ICardelh et al.l d 19891) with Ry = 3.1, is ap- 
plied to the synthetic spectral energy distributions calculated by 
the models to derive Lfir. The extinction is parameterized by 
the value of the color excess E(B-V). Since Lfir saturates for 
E(B-V) > 0.5, the value computed for E(B-V)= 1.0 should be 
taken as the maximum FIR emission expe cted in the starburst 
dMas-Hesse. Oti-Floranes & Cervinoll2008l) . We slightly mod- 
ified the extinction law, assuming that photons with a wave- 
length of 912 A < A < 1250 A suffered the same extinction 
^1250 A- I n tn i s wa y> we both a ) take into account photons with 
912 A < A < 1000 A, which would otherwise be ignored, and 
b) use a more realistic A,\ value, bearing in mind that the at- 
tenua tion law does not s eem to increase in this range (as Eq. 
(5) in lCardelh et al. ( 198 9) would imply), but to flatten, as con- 
cluded bv lMezger et al.l d 1982b . The effect of including the pho- 
tons with 912 A < A < 1000 A in the FIR calculation is that 
Lfir is increased by 5%. Similarly, the computed Lfir changes 
by up to 7% due to the modifications performed in the extinc- 
tion law in the range of912A</i< 1250 A, but this is only 
detectable for low E(B-V) values. A fraction 1 - / = 0.3 of ion- 
izing photons, irrespective of their energy, is assumed to be ab- 
sorbed bydust be fore they can ionize any atom, as derived by 
Mezger (1978) and Degi oia-Eastwood] d 1 9921) and recommended 
bv lBelfort et al.l dl987l) . The hypothetical presence of totally ob- 
scured stars was not considered. Since no assumption is made 
about the temperature of the dust, the models do not yield the 
infrared spectrum, but only the total energy absorbed by dust, 
which would be reemitted in the FIR range. 

The same prescriptions as those considered in CMHK02 
were taken with the SB99 models to model the FIR em ission 
of the burst. The extinction law by ICardelli et al.l (Qj89), with 
the minor changes explained above, was applied to the predicted 
SEDs, and 30% of the ionizing radiation was assumed to be con- 
verted into FIR emission. The estimation of Lfir using the SB99 
models is ~ 7% higher than predicted by CMHK02 for the same 



conditions, apparently due to the differences in the stellar atmo- 
spheres used, especially in the Lyman continuum . 



3.2. Ionizing power (N Lyc ) and emission lines 

As discussed above, we considered that a fraction 1 — / = 0.3 of 
the ionizing photons emitted by the stars is directly absorbed by 
dust, so that A^ yc represents only the number of photons which 
eventually could ionize the atoms in the surrounding gas. Since 
the intensities of the recom bination lines are directly propor- 
tional to A^Lyc (Osterbrock 1989), we wish to remark that our 
predictions are 30% lower than the luminosities computed by 
models which do not take this effect into account, which oth- 
erwise is realistic and should be included. Nevertheless, for the 
SC02 Pop. Ill predictions we have considered that no ionizing 
photons are destroyed, since at this extremely low metallicity no 
dust particles should as yet be present. 

The intensity of the emission lines has been derived from the 
predicted A^yc value, assuming Case B conditions (T e = 10 4 K, 
«e_= 500 cm 3 ) and the following relations dStorey & Hummeil 
fT995l) : 

L(Ha)= 1.36 x 10~ 12 N Lyc erg s" 1 
L(Ha)= 2.86 L(UB) 
L(Lya)= 8.7 L(Ha). 

We have computed the intensities of the emission lines for 
various E(B-V) values, but we have not considered the effects of 
radiation transfer on Lyor photons. It is well known that reso- 
nant scattering by neutral hydrogen significantly affects the Lya 
emission line profile and can lead to the total suppression of 
the line ( Verhamme et al. 2006), depending on the column den- 
sity and kinematics of the neutral gas an d the amount of dust 
dKunth et al.lll998l:lM~as-Hesse et al]|2003l) . The Lya luminosity 
we have computed therefore corresponds to the upper value ex- 
pected to escape from the star-forming region, corrected for the 
interstellar extinction in the same way as the Balmer lines, but 
not affected at all by neutral hydrogen scattering. In order to de- 
rive the star formation rate from the observed Lya luminosity 
this effect should be corrected before applying the calibration. 
It is important to note that the Lya photons escape fraction ca n 
cover a wide range, from to 1, as shown bv lAtek et al.l (|2009), 
though i n nearby starbur st galaxies the escape fraction is close 
to 10% (lOstlin et al.ll2008[). and around 5% in average in star- 
forming galaxies at z ~ 2 dHaves et1Ill2009l) . 



3.3. Continuum luminosity 

We have computed the luminosity evolution of the burst's con- 
tinuum emission at several wavelengths within the UV, optical 
and IR spectrum: 1500 A, 2000 A, 3500 A (U), 4400 A (B), 
5500 A (V) and 2.2 //m A (K). We included both the stellar and 
nebular continuum components and have considered the effect 
of interstellar extinction as parameterized by E(B-V). The most 
massive stars dominate the UV emission of star-forming regions 
and therefore drive the evolution of L1500 and L2ooo- On the other 
hand, the contribution of less massive stars can become impor- 
tant at longer wavelengths. In EB models, medium/low-mass 
stars accumulate due to their longer lifetimes and become the 
dominant contributors to L4400, L5500 and L22200 after ~ 100 Myr 
of evolution. 



6 



Otf-Floranes & Mas-Hesse: Calibration of star formation rate tracers 




. . . . Kennic:iitt(1998) 



Time (Myr) Time (Myr) 

Fig. 2. Evolution of L F j R predicted by CMHK02 for instantaneous (IB) models (left, normalized to SFS = 1 M ) and by SB99 for 
extended (EB) models (right, normalized to SFR = 1 M Q yr~'). Lfir has been computed for E(B-V)- 1 .0 and three metalli city values 
(from bottom to top: Z = 0.020-black, 0.008-blue, 0.001-red). We also plot on the right panel the Lfir calibration by iKennicuttl 
( 1998), adapted to the IMF used in this work. 
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Fig. 3. Evolution of A^yc predicted by CMHK02 for IB models (left, normalized to SFS - 1 M ) and for EB models (right, normal- 
ized to SFR = 1 M Q yr -1 ) as a function of metallicity (from bottom to to p: Z = 0.0 2 0-blac k, 0.008-blue, 0.004-red, 0.001-green), 
together with the / factor corrected and non-corrected A^yc values from IKennicuttl d 19981) . both adapted to the IMF used in this 
work. 



3.4. X-ray luminosity (L x ) 

The X-ray luminosity associated to a massive star formation 
episode is directly related to the amount of mechanical energy 
released into the interstellar medium (ISM) by stellar winds and 
supernova explosions. In order to calculate the mechanical en- 
ergy relased by supernova explosions, both the CMHK02 and 
SB99 models assume that each SN injects 10 51 erg into the ISM, 
but CMHK02 models subtract the energy emitted as X-ray ra- 
diation by SNR afterwards from the total ener gy of the super- 
nova. See Cervino, Mas -Hesse & Kunthl d2002l) for a more de- 
tailed description. Massive early-type stars, and especially Wolf- 
Rayet stars, lose a significant fraction of their mass in the form of 
strong winds. These winds end up interacting with the interstel- 
lar gas, which becomes extremely heated. As a result, soft X-ray 
emission is produced. Both models use the expressions of the 
gas terminal velocity from Leith erer. Robert & Drisser] (l992) 
together with their own prediction of mass loss to calculate the 
mechanical energy released by the winds. 

We have computed the evolution of the mechanical energy 
released by the burst per unit time, Ek- This energy is converted 
into X-ray luminosity assuming an efficiency factor e xe fr . As dis- 



cussed by iMas-Hesse. O tf-Floran es & Cervino! d2008l) the soft 
X-ray luminosity of star-forming galaxies can be reproduced by 
our synthesis models assuming e xe ff in the range of 1 - 10%. We 
assumed e xe ff= 0.05 for the calibrations presented in this work. 
The total soft X-ray luminosity is then obtained by adding the 
contribution by the SNR present in the star-forming region, still 
active as X-ray emitters. 



3.5. Radio luminosity (L md ) 

Both radio components, thermal and non-thermal, 
were calculated at 1.4 GHz. For the thermal 
emission L t hWe followed the presc riptions by 
iLequeux. Maucherat-Joubert. Deh arveng & Kunthl (|T98 1 ). 
This thermal component of the radio emission is directly 
proportional to the ionizing power of the starburst, with a slope 
art, = -0.1: 

La, = 1.57X 10- 26 (v(GHz))- 01 (T e /10 4 K) - 34 Af Lyc (erg s" 1 Hz" 1 ). 

We assumed an electronic temperature T e = 10 4 K and a 
frequency v = 1.4 GHz (i.e., 21.4 cm). 
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Fig. 4. Evolution of L t h predicted by CMHK02 for IB models (left, normalized to SFS = 1 M ) and for EB models (right, normalized 
to SFR = 1 M yr" 1 ) as a fu nction of metallicity (from bottom to top: Z = 0.020, 0.008, 0.004, 0.001). The dotted line shows the 
prediction by Condon] (1 19921) for EB models, adapted to the IMF used in this work. 
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Fig. 5. Evolution of L nth predicted by CMHK02 for IB models (left, normalized to SFS = 1 M Q ) and by SB99 for EB models (right, 
normalized to SFR = 1 M Q yr 1 ) fo r the metallici ty values Z = 0.020 (solid line), 0.008 (dashed line), 0.001 (point-dashed line). The 
dotted line shows the prediction by ICondonl (11992b for EB models, adapted to the IMF used in this work. 



The models by iMas-Hesse & Kunthl (Il99ll) underestimated 
the non-th e rmal emission associated to a starburst region. 
iMas-Hessd d 19921) showed that for the ratio L F]R /L t ac i to remain 
rather constant during the evolution of the burst, so that the ob- 
servational well-known correlation between radio and FIR emis- 
sion could be explained, the total non-thermal emission observed 
should be ~ 10 times larger than computed. It was argued that 
interactions between individual supernova remnants would yield 
a larger radio emission than observed in isolated SNR in our 
Galaxy, as the ones used f or the standard calibra tion. The same 
phenomenon was found bv lCondon & Yinl (|l990) when they ap- 
plied their Eq. (7) to our Galaxy, stating that the disagreement 
arises because the effect of accelerated electrons after ~ 2 x 10 4 
yr is ignored. They obtained their Eq. (8), which seems to over- 
come this problem. We used the expression 

L nth = 1.2 X 10 30 (v(GHz))-°- 9 v SN (erg s" 1 Hz" 1 ), 

where vsn is the supernova rate. This expression is similar 
to Eq. (8) from lCondon & Yinl dl99Qb . but assumes ff nt h = —0.9 
instead of a nt h = -0.8. For consistency, the factor used was 1.2 
instead of 1.3 to re produce the Galacti c val ues for L n ± at 40 8 
MHz and from iBerkhuii sen! (11984 and iTammannl (1 19821) , 
respectively. The L n ± value at 1.4 GHz obtained with this ex- 



pression is therefore 1 1 % lower than that from ICondon & Yinl 
(fl990l) . 

We have finally calibrated both SFS and SFR based on the 
total radio luminosity at 1 .4 GHz L rac j = L t h + L n th- 



4. Results and discussion 

In Tables[T|to[8]we present the calibration of the different estima- 
tors computed with our models. We list the values of the factor 
C A by which the different observables must be multiplied to ob- 
tain either SFS or S FR, i.e. SFS(L A ) = C A L A ,SFR(L A ) = C A L A 
As discussed above, M as-Hesse. Oti-Floranes & Cervifiol (f2008) 
showed that it is very important to take into account the evo- 
lutionary state of the burst when estimating its strength, espe- 
cially in the case of IB regimes, but also in some cases for 
EB episodes. For IB models we have calibrated the estima- 
tors for ages of 4, 5 and 6 Myr, which is the typical range of 
ages measured for star- formi ng galaxies (IMas-Hesse & K unth 
119991: iPellerin & Robert! [2007). The ionizing power decreases 
very rapidl y, so that the emission lines are barely detectable af- 
ter 7 Myr: IMas-Hesse & Kunthl (1 1999b predict EW(H/5) values 
below 2 A after this age for solar metallicity starbursts. 
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Fig. 6. Evolution (from top to bottom) of L1500, £2000 and L3500 predicted by CMHK02 for IB models (left, normalized to SFS = 1 
Mq) and by SB99 for EB models (right, normalized to SFR = 1 M yr" 1 ) both fo r Z = 0.020 (sol i d line s) and Z = 0.008 (dashed 
lines). Also in EB graph, values for Z4500 from the calibrations bv lKennicuttl d 1998b and lSalim et al.l (T2007), both adapted to the IMF 
used in this work. 



On the other hand, in an EB regime stars are replaced by new 
massive stars as they die, eventually reaching a steady state after 
a few tens of Myr. For this reason, we studied a wider range of 
ages: a non-evolved continuous episode with 10 Myr, a young 
burst after 30 Myr of evolution, and a star formation process al- 
ready in the steady state, at 250 Myr. Not all parameters remain 
constant when the steady state is reached, since medium/low- 
mass stars accumulate and contribute more and more to the op- 
tical - IR stellar continuum. As a result the calibration of some 
tracers is a function of the evolutionary state as well, though the 
effect is much weaker than for IB regimes. For A^yc and hence 
for L(Ha) and L(Lya), the asymptotic state is reached very soon, 
so that the same Ca value is displayed for the three points at 10 
Myr, 30 Myr and 250 Myr. 

For the parameters which are affected by interstellar extinc- 
tion we computed the calibration of the different estimators for 
different E(B-V) values within the range E(B-V)- 0.0-1.0 under 
the as sumption of the Galactic extinction law by Cardell iet alj 
(119891) . 

As mentioned above, our calibrations have been computed 
for a Salpeter initial mass function with limits at 2 and 120 M . 
To ease the comparison with other calibrations we list the cor- 
rection factors in Tables [TTIand[T2lbv which values from Tables 
HUH and|5][8] respectively, must be multiplied to convert our cal- 
ibrations to other usual IMF mass ranges of M — 0.1 - 100 M 
and M = 1 - 100 M Q . Similarly, in Tables [T3l and FBI we list 
the conversion factors which would be required for two limiting 
IMF slopes, with a = 1 and 3. These correction factors are the 
values L(2- 120, 2.35)/L(Mi ow -M h i g h, a), where M\ QVI and M high 
are the new mass limits, a is the new IMF slope and L the mag- 
nitude. For instance, the SFR(L ni a) calibration for a burst that has 
been creating stars with a constant rate for 30 Myr, and assum- 
ing an IMF with mass range M = 0.1 - 100 M Q and a = 2.35, is 
SFR(L mA ) = 3.51 4.3 x 10- 29 L rad = 1.5 X 10- 28 L rad (luminosities 
are measured in erg s _1 Hz T 1 and SFR in M yr _1 ). We calcu- 
lated the values of this correction factor using the SB99 models, 
since these models allow the IMF to be defined by the user. 

This correction between different IMFs depends on the age 
of the burst, on the parameter studied and on the star formation 
history assumed in the model and has two different origins: a) 
when considering different mass limits, both SFR and SFS mea- 
sure the mass converted into stars in a different mass range, and 



b) if Mhigh or a changes, so does the amount and eventual type 
of massive stars that contribute to the emission, which can be 
significantly affected. Emission by low-mass stars does not con- 
tribute dramatically to any of the computed parameters at the 
ages considered, and therefore any change in M\ ovl only exerts 
a variation of type a) in the calibration, as long as it does not 
increment unphysically. Also, since stars with masses above 100 
M Q have already died in IB models at the ages considered in the 
calibrations, only correction of type a) is needed (if a is not mod- 
ified), which can be performed an alytically by the calculation of 
the ra tio of the integrated IMFs dSchaererl 2003; Wi lkins et ail 
2008), and which does not depend on the evolutionary state. On 
the contrary, in EB models there are always stars more massive 
than 100 M G , hence a proper correction for both a) and b) us- 
ing SB99 as explained is needed. As can be checked in Tables 
[TTIand [T2"l and also b y a comparison of these values with those 
from Schaerer (2003), values of correction factors obtained an- 
alytically agree with those obtained for EB models using SB99 
typically within 5%, with the difference being higher only for 
A^Lyc (~ 10 - 15%) due to the fact that it is solely dominated 
by the most massive stars of the burst. Since SB99 models do 
not yield L so ftx, it was not possible to follow the same procedure 
in this case as for the other magnitudes, and therefore only the 
analytical value for L so ftx is given in Table [T2l 



The effect of the age of the burst in the correction of SFR 
estimators when different mass limits are assumed is very low 
(~ 8% at most) for the ages considered, hence we decided 
to leave only the asymptotical value at 250 Myr in Table [12] 
Similarly, the effect of E(B-V) value in the SFR(FIR) correc- 
tion factor is typically below 5%, therefore we opted for show- 
ing only the correction obtained for E(B-V) = 0.3, which lies 
somewhat in the middle of the range studied, where Lfir has not 
attained the saturation described in Sect. 13.11 Since L(Hor) and 
L(Lyof) are proportional to A^yo only the correction for the latter 
was included. We have applied these corr ections wh e n we com- 
pared our c a librati ons with the results of iKennicuttl (1 19981) and 
Sali m et al. (2007), which were comp uted for a Salp eter IMF, 
but with mass limits of 0.1 - 100 M .|Condon|(|l992) assumed 
mass limits of 5 - 100 M , but an IMF slope a = 2.50, there- 
fore L(2 - 120, 2.35) /L(5 - 100,2.50) was calculated explicitly 
in order to compare his radio calibrations with our predictions. 
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Our calibration of star formation estimators at different en- 
ergy ranges (from X-rays to radio) can also be used inversely 
(taking care of the different time constants and properties of the 
different tracers): for example, once the SFR or SFS is estimated 
from the UV or optical continuum luminosity, the expected 
emission at other ranges or emission lines can be directly com- 
puted, corrected for the effects of interstellar extinction in a con- 
sistent way. Nevertheless, the inverse use of the SFR, SFS cali- 
brators to derive the expected luminosity at other energy ranges 
has to be taken with caution. They are valid only for galaxies 
dominated by an ongoing massive star formation process, and 
so they would predict only the luminosity associated to the star 
formation episode. This could be just a fraction of the total lumi- 
nosity of the galaxy at other wavelength ranges dominated by an 
old underlying stellar population, like the near infrared. We have 
built a Web tool to compute the SFR or SFS for any given esti- 
mator and at different conditions, which is publicly accessible at 
http : //www . laeff . cab . inta-csic . es/research/sfr/ 

In the sections below we discuss the properties and applica- 
bility of each of the SFR and SFS estimators, and we compare 
our results with previous calibrations. 

4.1. Far infrared luminosity (Lfir) 

Evolution of Lfir is shown in Fig. [2] for both the CMHK02 
IB and SB99 EB models and different metallicities Z = 
0.020, 0.008, 0.001, computed for E(B-V)= 1.0. Since L FIR in 
star-forming regions is directly linked to the UV-optical emis- 
sion produced by the most massive stars of the bursts, it shows a 
peak in IB models at ~ 3 Myr, followed by a steep decline after- 
wards. On the other hand, when a constant SFR is assumed, Lfir 
starts to stabilize after the first 20 Myr of evolution, when an 
equilibrium between formation and destruction of massive stars 
is reached. Nevertheless, since medium/low-mass stars also con- 
tribute to Lfir, this parameter increases slowly but steadily as 
long as the star formation episode i s activ e. As discussed by 
Mas- Hesse. Oti-Floranes & Cervinol £2008), Lfir saturates for 
E(B-V) > 0.5 since the attenuation is already so high that most 
of the continuum emission is absorbed. We have calibrated Lfir 
as SFR and SFS estimator in Tables Q] to [8] for various E(B-V) 
values. If E(B-V)- 1.0 is assumed, Lfir becomes indeed a good 
estimate of the bolometric luminosity of the star formation pro- 
cess. 

The evolution of Lfir with time is sensitive to the metallicity, 
since low-metallicity stars evolve more slowly. As a result, Lfir 
is 1 5 - 30% larger for Z = 0.008 than for solar metallicity models 
at ages 4-6 Myr. In EB models, the effect is to accumulate 
more massive stars at a given time, so that Lfir becomes ~ 8 - 
1 1 % larger for Z = 0.008 once the burst is older than 30 Myr. 
On the other hand, during the initial phases of evolution, Lfir is 
essentially independent of the metallicity. 

Together with our predictions we show in Fig. [2] t he Lfir 
value derived from the expression by IKennicuttl d 1998b . which 
was IMF-corrected as explaine d above to be con sistent with the 
IMF we have used. The adapted K ennicutj dT998) value becomes 
log Lfir = 43.90, which agrees quite well with our predictions 
for solar metallicity at ages above ~ 100 Myr. On the other hand, 
the discrepancy becomes larger for less evolved bursts, differing 
by ~ 21% at 30 Myr, and by muc h more for the younger the 
episode. Although IKennicuttl d!998l) affirms that his calibration 
applies only to bursts younger than 100 Myr, it can be seen in 
Fig. |2] that it could be used accurately for solar metallicity ex- 
tended bursts at 50 - 250 Myr. 



4.2. Ionizing power (N Lyc ) and emission lines 

We show in Fig. [3] the evolution of A^ yc f° r IB an d EB models 
at different metallicities, as predicted by CMHK02 models. The 
emission of ionizing photons is driven exclusively by the most 
massive stars in the burst. Since these stars have short lifetimes, 
the ionizing power of instantaneous bursts decreases drastically 
with time after 3 Myr, when the most massive stars begin to ex- 
plode as supernovae. For the same reason, their number reaches 
an equilibrium value quite rapidly in models with a constant star 
formation rate, so that iV Lyc reaches an asymptotic value in just 
~ 8 Myr in EB models. This steady state of Nj_, yc lasts hundreds 
of Myr in EB models, but we show in the figure just the first 20 
Myr of evolution of the burst to remark its rapid stabilization. As 
we discussed, we considered that a fraction 1 — / = 0.3 of the 
ionizing photons are absorbed by dust, and here we show their 
number which can eventually ionize the atoms in the surround- 
ing gas. 

When considering Z < Zq, we find that the slower evolution 
of high-mass stars translates into a delay in the emission of ion- 
izing photons in IB models, as can be observed in Fig. [3] The 
effect is to increase A^yc significantly, by 0.4 dex at 4 Myr and 
0.7 dex at 6 Myr for Z = 0.001. In EB models the accumula- 
tion of massive stars originates A^yc values 0.2 dex higher for 
Z = 0.001. 

We include in Fig. |3] the N Lyc value by IKennicuttl d!998l) . 
adapted to our IMF limits. This A^yc value is larger than our 
prediction for solar meta llicity mod e ls by 50%, which is appar- 
ently due to the fact that IKennicuttl (1 19981) did not a pply any / 
correc tion. If we add 1 — / = 0.3 correction to the IKennicuttl 
( 1998) A^Lyc value, both predictions are within 6% in the asymp- 
totic phas e, as can be observed in Fig . [3] We insist that, as dis- 
cussed by Mas -Hesse & Kunthl (Il99ll) . it should be more realis- 
tic to assume that a fraction around 30% of ionizing photons are 
absorbed by dust, even in relatively dust-free environments, and 
that they therefore do not contribute to the ionization process. 

In Tables Q] to [8] we list the predicted A^yc values for IB and 
EB models at different evolutionary states and the predicted Ha 
and Lya calibrations for different values of interstellar extinc- 
tion, computed as discussed in Sect. 13.21 In Tables |9l and [TOl we 
list the corresponding A^yc predictions for starbursts dominated 
by Pop. Ill stars. As discussed above, 1 — / = 0.0 has been as- 
sumed in this case. 

We wish to remark that there may be significant differences 
between the ex tinction affect i ng the nebular lines and the stel- 
lar continuum. ICalzetti et al.l d2000l) estimated that on average 
the color excess (E(B-V)) derived from the stellar continuum 
was a factor of 0.4 lower than the value derived from the emis- 
sion lines. The analysis of NGC 4214, a well-resolved starburst 
galaxy, by Maiz-Apella niz etaf] ([1998) showed that the dust ap- 
peared concentrated at the boundaries of the ionized region in 
its central star-forming region and affected mainly the nebular 
emission lines, while the stellar continuum itself was located in 
a region relatively free of dust and gas. But on the other hand, 
no spatial decoupling between massive stars, ionized gas and 
dust was found on other star-forming regions of the same galaxy, 
which indicates that the specific geometry of each star-forming 
knot driv es the decouplin g of continuum and emission line ex- 
tinctions. [ErbetjO] (|2006|) concluded from the analysis of a large 
sample (114 objects) of star-forming galaxies at z ~ 2 that the 
average SFR values derived from both L(Hc) and Ljjy, apply- 
ing the same E(B-V) for the extinction correction, were com- 
pletely consistent (indeed, almost identical). Nevertheless, the 
scatter was significant (0.3 dex), which again indicates that the 
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Fig. 9. Evolution of L22200 predicted by SB99 for EB models as 
in Fig. |g] 



geometry of individual objects plays an important role. This ef- 
fect has to be taken into account if the SFR is derived simulta- 
neously from tracers based on both the emission lines intensity 
and continuum level. Finally, we want to stress that the strength 
of the extinction assumed in our computations re sults from the 
comb ination of the listed E(B-V) values and the ICardelli et alj 
(1989) law. To get the same extinction with different laws, cor- 
respondingly different E(B-V) values have to be used (for ex- 
ample, the same extinction i s obtained at Ha with E(B-V)= 0.3 
and the ICardelli et all dl989l) law as with E(B-V)= 0.23 and the 
ICalzetti et al.l (120001) parameterization). 

4.3. Radio emission 

The evolution of L t h and L nt h is presented in Figs. Hand |3J re- 
spectively. Since the thermal radio component is assumed to be 
proportional to A^yc, its evolution directly follows that of A^yc, 
as discussed above. But we still show the plots in Fig.|4]to allow 
for a direct comparison between the intensities of the thermal 
and non-thermal components at different evolutionary states of 
the burst. We have also p lotted on Figs. [4] and [5] the values pre- 
dicted bv lCondonl (Q~992) for EB models after they were adapted 
to our IMF. 

When constant star formation is assumed, L nt h is observed 
to increase two orders of magnitude within ~ 35 Myr after the 
SN start to explode, reaching then a steady state with almost no 
sensitivity to the metal content of the star-forming region. The 
predictions for IB models, on the other hand, are more sensitive 
to metallicity, especially at ages between 3 and 12 Myr, when the 
most massive stars explode (their lifetime and evolution depend 
on metallicity). 

In Tables[T]to[8]we list the calibration of Ln&, which includes 
both the thermal and non-thermal components at any time. It can 
be seen in the figures that the total radio emission at 1 .4 GHz is 
dominated by the non-thermal component as soon as supernova 
explosions take place in the star-forming region. Only during the 
first 2-3 Myr in IB bursts can the radio luminosity be dominated 
by thermal emission. 

Although the calibration derived by ICondonl d 19921) is con- 
sistent with our predictions, Fig. |4] shows that the value is 
12% higher than predict ed by our mode ls. However, when using 
the L(Ho-) value from IKennicuttl d 19831) as ICondonl d 19921) did, 
with the prescriptions we have explained in Sect. 13.21 it can be 
checked that th is difference falls to 5%. But we wish to note that 
ICondonl dl992h ignores the dust absorption of ionizing photons, 



but assumes a ~ 30% lower ionizing power of the burst than 
predicted by the SB99 models. Both effects seem to cancel each 
other, which renders both L t h predictions consistent. 

4.4. Ultraviolet and U-band continuum 

UV emission is dominated by the most massive stars in the burst 
(via stellar and nebular continuum radiation), and therefore con- 
stitutes a direct tracer of the presence of young, massive stars, 
since it is free of contamination by the continuum emission 
of low-mass stars. Nevertheless, interstellar extinction is very 
strong in this energy range, and therefore its correction plays a 
key role in the use of the UV continuum luminosity as a tracer 
of star formation. In Fig. [6] we show the evolution of the lumi- 
nosity at 1500, 2000 and 3500 A for IB models as predicted 
by CMHK02, and for EB models by SB99, for Z = 0.020 and 
Z = 0.008. No extinction has been applied in these figures. 

In IB models, the UV continuum follows the well-known 
trend, with a peak of the emission at ~ 3 Myr, and a severe de- 
crease when the most massive stars begin to explode. The slower 
evolution at low metallicity is also evident in the figures. In the 
case of EB models, the UV continuum asymptotically stabilizes 
after the first ~ 50 Myr of evolution. The stabilization is slower 
for L3500, since it is contributed by stars of a lower mass than 

Ll500- 

We show also in Fi g.|6]fhe L1500 values bv lKennicuttl (119981) 
and lSalim et al.l (120071) . after they are adapted to our IMF mass 
limits. Both calibrations ar e close to our pre dictions, with a per- 
fect match with the one bv lKennicuttl d!998l) at ages around 100 
Myr fo r solar metallicity models. The predictions bv lSalim et al.l 
(20071) are so mewhat ov e restim ated with respect to both the 
predictio ns by IKennicuttl d l998) and by us. But these authors 
used the Bruzual & Chariot (2003) synthesis models in order 
to fit the photometric data points of their sample, assuming an 
exponentially-declining continuous star formation history, over 
which random-in-time starbursts were superimposed. The star 
formation histories are so different between these models that a 
direct comparison is not straightforward. 

IKennicuttl (1 19981) argued that his calibration is applicable to 
bursts which have been forming stars for 100 Myr or longer. Our 
results show that these calibrations can indeed be used accurately 
for ages above ~ 30 Myr, since the stronger source of error will 
be the correction of the interstellar extinction. 

In Tables |9land[T0l we list the corresponding L1500 predictions 
for starbursts dominated by Pop. Ill stars. As discussed above, no 
extinction has been assumed in this case, since no dust particles 
are expected to be present. 

4.5. B-, V- and K-band continuum 

In the optical and IR bands the contribution to the continuum 
emission by low-mass stars becomes very important as the star- 
burst evolves. In Fig. [7] we show the evolution of both L4400 and 
L5500 for different metallicities. In IB models the formation of 
Wolf-Rayet stars induces some bumps between 3 and 6 Myr, 
due to their short lifetime and their relatively low number. These 
bumps are indeed a stochastic effect and should not be consid- 
ered as real features. The effect of low-mass stars accumula- 
tion is clearly seen in the curves for EB models. The slope is 
higher than for the continuum at shorter wavelengths and, in- 
deed, L4400 and L5500 do not stabilize even after 250 Myr of 
evolution. Moreover, the continuum in the optical range can be 
strongly contaminated by an underlying older stellar population, 
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Fig. 7. Evolution (from top to bottom) of L4400 and L5500 predicted by CMHK02 for IB models (left, normalized to SFS = 1 M ) 
and SB99 for EB models (right, normalized to SFR = 1 M yr _1 ) both for Z = 0.020 (solid lines) and Z = 0.008 (dashed lines). 
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Fig. 8. Evolution of L22200 predicted by CMHK02 for IB models (left, normalized to SFS = 1 M ) and SB99 for EB models (right, 
normalized to SFR = 1 M yr 1 ). From top to bottom: Z = 0.020-black, 0.008-blue, 0.001-red. 
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Fig. 10. Evolution of E K predicted by CMHK02 for IB models (left, normalized to SFS = 1 M ) and by SB99 for EB models (right, 
normalized to SFR = 1 M yr _1 ) as a function of metallicity. From top to bottom: Z = 0.020-black, 0.008-blue, 0.001-red. 



not related to the present episode of star formation, which makes 
the estimators based on L4400 and L5500 very uncertain. 

The evolution of L22200 is strongly affected by the formation 
and evolution of red supergiant stars (RSG), which are very de- 
pendent on metallicity. We plot in Fig. [8] the evolution of L22200 
as predicted by CMHK02 for IB models and by SB99 for EB 
models for three values of the metallicity. The formation of RSG 



is favoured in high metallicity environments, dominating com- 
pletely L22200 in Z = 0.020 IB models in the period of 5 - 15 
Myr. On the other ha nd, their contribution for Z = 0.001 be- 
comes negligible (see ICervino &Mas-HesseNl994 for a more 
detailed discussion). 

In extended episodes the population of RSG stars tends to 
stabilize at high metallicities. Nevertheless, the accumulation of 
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low-mass stars, whose contribution to the K-band is important, 
drives a continuous increase of L22200, which is specially evident 
for low metallicities where the contribution of RSG is negligible. 
In Fig. [9] we zoom on the first 20 Myr of evolution, to show the 
effect of metallicity on L22200 early evolution. 

4.6. X-ray luminosity (L x ) 

At the early stages of the burst, the injection of mechanical en- 
ergy into the ISM is dominated by stellar winds, which tend to 
be stronger as the most massive stars evolve. Around 3 Myr after 
the onset of the burst, massive stars end their lifetimes, and this 
component begins to decrease monotonically. Then the first SN 
start to explode, becoming the dominant sources of mechanical 
energy release in IB models after 3-4 Myr. In EB models both 
components tend to be similar once the stellar population has at- 
tained a steady state, around 30 - 40 Myr after the onset of the 
burst. This can be seen in Fig. [10] where we have plotted the evo- 
lution of E K as predicted by SB99 and CMHK02 for EB and IB 
models respectively, and for different metallicities. The strongest 
release of mechanical energy in IB models occurs at ~ 4 Myr, 
when the first SN appear and add their contribution to the winds 
dominated by the Wolf-Rayet stars. The release of mechanical 
energy stabilizes in IB models after ~ 6 Myr of evolution, since 
the supernova rate decreases very slowy for the n ext 20 Myr for 
a Salpeter IMF (see lCervifio & Mas-Hessd (11994 ). In EB mod- 
els the evolution of Ek stabilizes completely 30 - 40 Myr after 
the onset of the burst, depending on metallicity, when both the 
population of the most massive stars, responsible for the winds 
and SN, and medium-mass stars (M > 6 - 8 M ), originator of 
SN as well, have reached an equilibrium. 

Ek is very sensitive to metallicity as long as it is dominated 
by stellar winds, since the mass-loss rate varies largely with the 
metal abundance of stellar atmospheres. This is especially evi- 
dent in the first 3 Myr of evolution of IB models, as shown in 
Fig. [10] but the effect can also be seen on the lowest metallicity 
EB models. 

As discussed above we have assumed e xe s= 0.05 to cal- 
ibrate Lsoftx, i-e., we consider that 5% of the mechanical 
power that is released is reconverted into soft X-ray emis- 
sion. We want to stress that while e xe s seems to be in the 
range e xe ff= 0.01 - 0.1 for star-forming galaxies, it could 
attain higher or lower values in individual objects, depend- 
ing on the details of the interaction between the accelerated 
winds and the surrounding interstellar medium. In Table [JJ 
we list only the S FS (£ so ftx) calibration for an IB model at 
5 My r. As discussed by iMas-HesseT Otf-Floranes & Cervinol 
d2008l) there are large uncertainties in L son x at these ages, there- 
fore the value listed should be considered as the best aver- 
age ^jOTjhe_4_j^_6_^My^_ger^ SFR calibration 
by iMas-Hesse. Otf-Floranes & Cervinol (120081) did not provide 
a value for a burst age of 250 Myr, but as we have explained, 
a steady state is expected at 30 - 40 Myr in EB models, so we 
assume the same estimator value for 250 Myr as for 30 Myr, as 
listed in Table [5] 

5. Conclusions 

We have computed a self-consistent calibration of the most usual 
star formation rate tracers using evolutionary synthesis models 
optimized for both short-lived and extended star formation pro- 
cesses. We have taken into account the effects of the evolutionary 
state, the star formation regime and the effect of interstellar ex- 
tinction on the different estimators. We have also included the 



prediction for population III starbursts. The main results can be 
summarized as follows: 

- The intensity of star formation episodes should be measured 
as a star formation rate (SFR), in units of M yr _1 , when star 
formation proceeds at nearly constant rates for long periods 
of time (tens to hundreds of Myr), and as a star formation 
strength (SFS), in units of M , when the episode is nearly 
instantaneous or short-lived (few Myr). Deriving an SFR 
value for an instantaneous burst would not have any phys- 
ical meaning. Moreover, the evolution is so rapid in these 
conditions that the results would be misleading and even er- 
roneous. 

- We have provided calibrations for SFS estimators at the ages 
of 4 - 6 Myr. The ionizing power of these bursts decays so 
rapidly that they would show just very weak emission lines 
at older ages, and would not be easily recognized as massive 
starbursts. 

- The evolutionary state should also be taken into account 
even for extended episodes of star formation, whenever pos- 
sible. In general most of the SFR tracers have stabilized after 
100 Myr of evolution at a constant rate. 

- Some of the most used SFR and SFS tracers are strongly af- 
fected by interstellar extinction. We have provided their cali- 
bration as a function of E(B-V) assuming the Galactic extinc- 
tion law bv lCardelli et alj d!989l) . to allow for a more precise 
result if the extinction can be somehow estimated. SFR(FIR) 
has also been computed as a function of the extinction. For 
E(B-V) values above 0.5 Lfir can be considered as a good 
estimate of the bolometric luminosity of the burst. 

- SFR and SFS tracers allow to derive the amount of gas trans- 
formed into stars, but they are calibrated for a given initial 
mass function. The choice of different usual IMF mass lim- 
its within the range 0.1 - 120 M , and IMF slopes a = 1, 3 
can lead to variations in the derived rates of up to a factor 
4. A comparison of results by different authors is therefore 
meaningful only after a normalization to a common IMF. 

- Our calibration of star formation estimators at different en- 
ergy ranges (from X-rays to radio) can also be used inversely, 
(taking care of the different time constants and properties of 
the different tracers): for example, once the SFR or SFS is 
estimated from the UV or optical continuum luminosity, the 
expected emission at other ranges or emission lines can be 
directly computed, corrected for the effects of interstellar ex- 
tinction in a consistent way. Nevertheless, the inverse use of 
the SFR and SFS calibrators to derive the expected luminos- 
ity at other energy ranges has to be taken with caution, since 
they would predict only the luminosity associated to the star 
formation episode, which could be just a fraction of the total 
luminosity of the galaxy at certain wavelength ranges domi- 
nated by an old underlying stellar population. 
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Table 1. SFS estimators which do not depend on E(B-V). 



Magnitude 


IB (4 Myr) 


IB (5 Myr) 


IB (6 Myr) 


N L yc (*"' Mo" 1 ) 

L^Cergs- 1 Hz-'Me-') 
^ so ftx (erg s" 1 Mq" 1 ) 


6.7 x 1(T 47 
7.9 x 10~ 22 


1.0 x 1(T 46 
7.9 x 10- 22 
2 x 10~ 34 


4.5 x 10~ 46 
9.5 x 10- 22 
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Table 2. SFS estimators which do depend on E(B-V) for an IB model at age=4 Myr. 
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E(B-V) 


L(Lya) 




L(Haf 




7 b 

-M500 




^2000 b 




I b 

i-3500 




7 b 

L4400 




7 b 

i-5500 




^22200 b 




^FIR a 




0.0 


5.7 x 10" 




4.9 x 10 




3.3 x 


io- 




6.7 x 


io- 




2.8 x 


io- 




4.2 x 


io- 




8.3 x 


io- 




7.8 x 


io- 




- 




0.1 


1.5 x 10 


-35 


6.2 x 10 


-35 


7.0 x 


io- 


34 


1.5 x 


io- 


-33 


4.4 x 


io- 


-33 


6.1 x 


io- 


-33 


1.1 X 


io- 


-32 


8.1 x 


io- 


31 


3.0 x 10 


-37 


0.2 


3.8 x 10~ 


35 


7.8 x 10 


-35 


1.5 x 


io- 


33 


3.4 x 


io- 


-33 


7.0 x 


io- 


-33 


8.9 x 


io- 


-33 


1.5 x 


io- 


-32 


8.3 x 


io- 


31 


2.1 x 10 


-37 


0.3 


9.9 x 10 


35 


9.9 x 10 


-35 


3.2 x 


io- 


33 


7.7 x 


io- 


-33 


1.1 X 


io- 


-32 


1.3 x 


io- 


-32 


2.0 x 


io- 


-32 


8.6 x 


io- 


31 


1.9 X 10 


-37 


0.4 


2.6 x 10- 


-34 
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-34 


6.9 x 


io- 


33 


1.7 x 


io- 


-32 


1.7 x 
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-32 


1.9 x 
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2.6 x 
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8.9 x 
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-37 


0.5 


6.6 x 10- 


34 


1.6 x 10 


-34 


1.5 x 


io- 


32 


3.9 x 


io- 


-32 


2.7 x 


io- 


-32 


2.8 x 


io- 


-32 


3.5 x 


io- 


-32 


9.2 x 


io- 


31 


1.7 x 10 
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8.1 x 
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-32 


1.0 x 


io- 


30 


1.7 x 10 


-37 
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3.0 x 10- 


-32 


4.0 x 10 


-34 


3.1 x 


io- 


31 


1.0 x 


io- 


-30 


1.7 x 


io- 


-31 


1.3 x 


io- 


-31 


1.1 X 


io- 


-31 


1.0 x 


io- 


30 


1.7 x 10 


-37 
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7.8 x 10- 


-32 


5.1 x 10" 


-34 


6.6 x 


io- 


31 


2.3 x 


io- 


-30 


2.7 x 


io- 


-31 


1.8 x 


io- 


-31 


1.4 x 


io- 


-31 


1.1 X 


io- 


30 


1.7 x 10 


-37 



a Measured in erg s 1 M Q 1 . 
b Measured in erg s -1 A -1 M ( 



Table 3. SFS estimators which do depend on E(B-V) for an IB model at age=5 Myr. 
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-33 


3.4 x 10- 
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4.7 x 


io- 


-32 


4.0 x 


io- 


-32 


3.7 x 


io- 


-32 


4.0 x 10- 


-31 


2.6 x 10 


-37 


0.6 


2.6 x 10- 


-33 


3.0 x 10 


-34 


4.8 x 


io- 


32 


1.4 x 10- 


31 


7.5 x 


io- 


-32 


5.9 x 


io- 


-32 


5.0 x 


io- 


-32 


4.1 x 10- 


-31 


2.6 x 10- 


-37 


0.7 


6.7 x 10- 


-33 


3.8 x 10 


-34 


1.0 x 


io- 


31 


3.1 x 10- 


31 


1.2 x 


io- 


-31 


8.6 x 


io- 


-32 


6.6 x 


io- 


-32 


4.2 x 10- 


-31 


2.6 x 10- 


-37 


0.8 


1.7 x 10- 


-32 


4.8 x 10 


-34 


2.2 x 


io- 


31 


7.0 x 10- 


-31 


1.9 x 


io- 


-31 


1.3 x 


io- 


-31 


8.8 x 


io- 


-32 


4.4 x 10- 


-31 


2.6 x 10 


-37 


0.9 


4.5 x 10 


-32 


6.1 x 10 


-34 


4.7 x 


io- 


31 


1.6 x 10- 


-30 


2.9 x 


io- 


-31 


1.8 x 


io- 


-31 


1.2 x 


io- 


-31 


4.5 x 10- 


-31 


2.6 x 10- 


-37 


1.0 


1.2 x 10" 


-31 


7.7 x 10- 


-34 


1.0 x 


io- 


30 


3.5 x 10- 


-30 


4.6 x 


io- 


-31 


2.7 x 


io- 


-31 


1.6 x 


io- 


-31 


4.7 x 10- 


-31 


2.5 x 10- 


-37 



a Measured in erg s 1 M Q 1 . 
b Measured in erg s~' A~' M _1 . 



Table 4. SFS estimators which do depend on E(B-V) for an IB model at age=6 Myr. 



E(B-V) 


L(Lya) 


a 


L(Ha)' 




7 b 

-M500 




^2000 b 




I b 

i-3500 




I b 

i-4400 




I b 

i-5500 




i-22200 b 




^FIR 




0.0 


3.8 x 10- 


-35 


3.3 x 10- 


-34 


6.9 x 


io- 


34 


1.5 x 


io- 


-33 


7.0 x 


10" 


-33 


5.8 x 


10" 


-33 


8.6 x 


10" 


-33 


2.5 x 


10" 


31 






0.1 


9.9 x 10- 


-35 


4.2 x 10 


-34 


1.5 x 


io- 


33 


3.3 x 


io- 


-33 


1.1 X 


io- 


-32 


8.4 x 


10" 


-33 


1.1 X 


10" 


-32 


2.6 x 


10" 


31 


6.7 x 10 


-37 


0.2 


2.6 x 10- 


-34 


5.3 x 10 


-34 


3.2 x 


io- 


33 


7.4 x 


io- 


-33 


1.7 x 


io- 


-32 


1.2 x 


10" 


-32 


1.5 x 


10" 


-32 


2.7 x 


10" 


31 


4.6 x 10 


-37 


0.3 


6.7 x 10- 


-34 


6.7 x 10 


-34 


6.7 x 


10" 


33 


1.7 x 


10" 


-32 


2.7 x 


10" 


-32 


1.8 x 


10" 


-32 


2.0 x 


10" 


-32 


2.7 x 


10" 


31 


4.0 x 10- 


-37 


0.4 


1.7 x 10- 


-33 


8.4 x 10 


-34 


1.4 x 


io- 


32 


3.8 x 


io- 


-32 


4.3 x 


10" 


-32 


2.6 x 


10" 


-32 


2.7 x 


10" 


-32 


2.8 x 


io- 


31 


3.8 x 10- 


-37 


0.5 


4.5 x 10- 


-33 


1.1 x 10 


-33 


3.1 x 


io- 


32 


8.5 x 


io- 


-32 


6.8 x 


10" 


-32 


3.8 x 


10" 


-32 


3.6 x 


10" 


-32 


2.9 x 


10" 


31 


3.7 x 10- 


-37 


0.6 


1.2 x 10- 


-32 


1.3 x 10 


-33 


6.6 x 


io- 


32 


1.9 x 


io- 


-31 


1.1 X 


10" 


-31 


5.6 x 


10" 


-32 


4.8 x 


10" 


-32 


3.0 x 


io- 


31 


3.6 x 10 


-37 


0.7 


3.0 x 10- 


-32 


1.7 x 10 


-33 


1.4 x 


io- 


31 


4.3 x 


io- 


-31 


1.7 x 


10" 


-31 


8.2 x 


10" 


-32 


6.3 x 


10" 


-32 


3.1 x 


io- 


31 


3.5 x 10" 


-37 


0.8 


7.8 x 10- 


-32 


2.1 x 10 


-33 


3.0 x 


io- 


31 


9.7 x 


io- 


31 


2.6 x 


10" 


-31 


1.2 x 


10" 


-31 


8.4 x 


10" 


-32 


3.2 x 


10" 


31 


3.5 x 10- 


-37 


0.9 


2.0 x 10- 


-31 


2.7 x 10 


-33 


6.5 x 


io- 


31 


2.2 x 


io- 


-30 


4.2 x 


10" 


-31 


1.7 x 


10" 


-31 


1.1 X 


10" 


-31 


3.3 x 


io- 


31 


3.5 x 10" 


-37 


1.0 


5.3 x 10- 


-31 


3.4 x 10 


-33 


1.4 x 


io- 


30 


4.9 x 


io- 


-30 


6.6 x 


10" 


-31 


2.5 x 


10" 


-31 


1.5 x 


10" 


-31 


3.4 x 


io- 


31 


3.5 x 10- 


-37 



a Measured in erg s 1 M Q 1 . 
b Measured in erg s~' A~' M ~'. 
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Table 5. SFR estimators which do not depend on E(B-V). 



Magnitude 


EB (lOMyr) 


EB (30 Myr) 


EB (250 Myr) 


Atyc^'CMoyr- 1 )-') 


4.1 x 10~ 54 


4.1 x 10~ 54 


4.1 x 10~ 54 


L rad (ergs-' Hz- 1 (M G yr" 1 )" 1 ) 


9.3 x lfr 29 


4.3 x 10- 29 


3.8 x 10- 29 


^softx (erg s- 1 (Moyr-'r 1 ) 


8 x 10- 41 


2x 10- 41 


2x 10- 41 



Otf-Floranes & Mas-Hesse: Calibration of star formation rate tracers 
Table 6. SFR estimators which do depend on E(B-V) for an EB model at age=10 Myr. 
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E(B-V) 


L(Lya) 


a 


L(Hq-) 1 




M500 




^2000 b 




I b 

i-3500 




I b 

L4400 




I b 

i-5500 




^22200 b 




f a 




0.0 


3.5 x 10 


43 


3.0 X 10 


-42 


4.1 x 


10" 


-4 1 


8.4 x 


10" 


41 


3.4 x 


10" 


-40 


6.5 x 


10" 


-40 


1.1 X 


10" 


-39 


2.3 x 


10" 


38 








0.1 


9.0 x 10- 


43 


3.8 x 10 


-42 


8.8 x 


10" 


-4 1 


1.9 x 


10" 


40 


5.3 x 


10" 


-40 


9.5 x 


10" 


-40 


1.5 x 


10" 


-39 


2.4 x 


10" 


38 


3.3 x 


10" 


-44 


0.2 


2.3 x 10" 


42 


4.8 x 10 


-42 


1.9 x 


10" 


-40 


4.3 x 


10" 


40 


8.4 x 


10" 


-40 


1.4 x 


10" 


-39 


1.9 x 


10" 


-39 


2.5 x 


10" 


38 


2.4 x 


io- 


-44 


0.3 


6.1 x 10~ 


42 


6.1 x 10 


-42 


4.0 x 


10" 


-40 


9.6 x 


10" 


40 


1.3 x 


10" 


-39 


2.0 x 


10" 


-39 


2.6 x 


10" 


-39 


2.6 x 


10" 


38 


2.2 x 


io- 


-44 


0.4 


1.6 x 10- 


41 


7.7 x 10 


-42 


8.6 x 


10" 


-40 


2.2 x 


10" 


-39 


2.1 x 


10" 


-39 


2.9 x 


10" 


-39 


3.4 x 


10" 


-39 


2.6 x 


10" 


38 


2.0 x 


io- 


-44 


0.5 


4.1 x 10- 


41 


9.7 x 10 


-42 


1.8 x 


10" 


-39 


4.9 x 


10" 


-39 


3.3 x 


10" 


-39 


4.3 x 


10" 


-39 


4.6 x 


10" 


-39 


2.7 x 


10" 


38 


2.0 x 


io- 


-44 


0.6 


1.1 x 10" 


40 


1.2 x 10" 


-41 


4.0 x 


10" 


-39 


1.1 X 


10" 


-38 


5.2 x 


10" 


-39 


6.3 x 


10" 


-39 


6.1 x 


10" 


-39 


2.8 x 


10" 


38 


2.0 x 


io- 


-44 


0.7 


2.7 x 10" 


40 


1.5 x 10 


-41 


8.5 x 


10" 


-39 


2.5 x 


10" 


-38 


8.2 x 


10" 


-39 


9.1 x 


10" 


-39 


8.1 x 


10" 


-39 


2.9 x 


10" 


38 


1.9 x 


io- 


-44 


0.8 


7.1 x 10" 


4(1 


2.0 x 10 


-41 


1.8 x 


10" 


-38 


5.6 x 


10" 


-38 


1.3 x 


10" 


-38 


1.3 x 


10" 


-38 


1.1 X 


10" 


-38 


3.0 x 


10" 


38 


1.9 x 


io- 


-44 


0.9 


1.8 x 10" 


-39 


2.5 x 10 


-41 


3.9 x 


10" 


-38 


1.3 x 


10" 


-37 


2.0 x 


10" 


-38 


1.9 x 


10" 


-38 


1.4 x 


10" 


-38 


3.1 x 


10" 


38 


1.9 x 


io- 


-44 


1.0 


4.8 x 10" 


-39 


3.1 x 10" 


-41 


8.3 x 


10" 


-38 


2.8 x 


10" 


-37 


3.2 x 


10" 


-38 


2.8 x 


10" 


-38 


1.9 x 


10" 


-38 


3.2 x 


10" 


38 


1.9 x 


io- 


-44 



a Measured in erg s 1 (M Q yr 1 ) 1 . 
b Measured in erg s -1 A~' (M Q yr -1 )" 1 . 



Table 7. SFR estimators which do depend on E(B-V) for an EB model at age=30 Myr. 



E(B-V) 


L(Lya) 


a 


L(Ha)' 




I b 

-M500 




^2000 b 




/ b 

i-3500 




/ b 

i-4400 




/ b 

i-5500 




i-22200 b 




1 ' d 
^FIR 




0.0 


3.5 x 10- 


43 


3.0 x 10 


-42 


3.3 x 


io- 


-41 


6.6 x 10- 


-41 


2.6 x 


10" 


-40 


4.2 x 


10" 


-40 


7.2 x 


io- 


-40 


1.2 x 10- 


38 








0.1 


9.0 x 10- 


43 


3.8 x 10 


-42 


7.1 x 


io- 


-41 


1.5 x 10- 


-40 


4.1 x 


10" 


-40 


6.1 x 


io- 


-40 


9.6 x 


io- 


-40 


1.3 x 10- 


38 


2.7 x 


io- 


-44 


0.2 


2.3 x 10- 


42 


4.8 x 10- 


-42 


1.5 x 


io- 


-40 


3.3 x 10" 


-40 


6.5 x 


10" 


-40 


8.9 x 


io- 


-40 


1.3 x 


io- 


-39 


1.3 x 10- 


38 


2.0 x 


io- 


-44 


0.3 


6.1 x 10- 


42 


6.1 x 10 


-42 


3.3 x 


io- 


-40 


7.5 x 10- 


-40 


1.0 x 


10" 


-39 


1.3 x 


io- 


-39 


1.7 x 


io- 


-39 


1.4 x 10- 


38 


1.8 x 


10" 


-44 


0.4 


1.6 x 10- 


41 


7.7 x 10 


-42 


7.0 x 


io- 


-40 


1.7 x 10- 


-39 


1.6 x 


10" 


-39 


1.9 x 


io- 


-39 


2.2 x 


io- 


-39 


1.4 x 10- 


38 


1.7 x 


io- 


-44 


0.5 


4.1 x 10- 


41 


9.7 x 10 


-42 


1.5 x 


10 


-39 


3.8 x 10- 


-39 


2.5 x 


10" 


-39 


2.8 x 


io- 


-39 


3.0 x 


io- 


-39 


1.5 x 10- 


38 


1.6 x 


io- 


-44 


0.6 


1.1 x 10- 


4(1 


1.2 x 10 


-41 


3.2 x 


io- 


-39 


8.6 x 10- 


-39 


4.0 x 


10" 


-39 


4.1 x 


io- 


-39 


4.0 x 


io- 


-39 


1.5 x 10- 


38 


1.6 x 


io- 


-44 


0.7 


2.7 x 10- 


4(1 


1.5 x 10 


-41 


6.8 x 


io- 


-39 


1.9 x 10- 


-38 


6.3 x 


10" 


-39 


5.9 x 


io- 


-39 


5.3 x 


io- 


-39 


1.6 x 10 


38 


1.6 x 


io- 


-44 


0.8 


7.1 x 10- 


40 


1.9 x 10 


-41 


1.5 x 


io- 


-38 


4.4 x 10- 


-38 


9.9 x 


10" 


-39 


8.6 x 


io- 


-39 


7.0 x 


io- 


-39 


1.6 x 10 


38 


1.6 x 


io- 


-44 


0.9 


1.8 x 10- 


-39 


2.5 x 10 


-41 


3.1 x 


io- 


-38 


9.8 x 10" 


-38 


1.6 x 


10" 


-38 


1.3 x 


io- 


-38 


9.4 x 


io- 


-39 


1.7 x 10- 


38 


1.6 x 


io- 


-44 


1.0 


4.8 x 10- 


-39 


3.1 x 10 


-41 


6.7 x 


io- 


-38 


2.2 x 10" 


-37 


2.4 x 


10" 


-38 


1.8 x 


io- 


-38 


1.2 x 


io- 


-38 


1.7 x 10" 


38 


1.5 x 


io- 


-44 



a Measured in erg s 1 (M Q yr 1 ) 1 . 
b Measured in erg s~' A~' (M Q yr -1 ) 



Table 8. SFR estimators which do depend on E(B-V) for an EB model at age=250 Myr. 



E(B-V) 


L(Lya) 


a 


L(Ha)' 




/ b 

-M500 




^2000 b 




I b 

i-3500 




I b 

i-4400 




I b 

i-5500 




^-22200 b 




7 a 
^FIR 




0.0 


3.5 x 10- 


-43 


3.0 x 10- 


-42 


2.7 x 


io- 


-41 


5.1 x 


io- 


-41 


1.8 x 


io- 


-40 


2.2 x 


io- 


-40 


3.6 x 


io- 


-40 


6.7 x 


io- 


39 








0.1 


9.0 x 10- 


-43 


3.8 x 10 


-42 


5.8 x 


io- 


-41 


1.2 x 


io- 


-40 


2.9 x 


io- 


-40 


3.2 x 


io- 


-40 


4.8 x 


io- 


-40 


6.9 x 


io- 


39 


2.3 x 


io- 


-44 


0.2 


2.3 x 10- 


42 


4.8 x 10 


-42 


1.2 x 


io- 


-40 


2.6 x 


io- 


-40 


4.5 x 


io- 


-40 


4.7 x 


io- 


-40 


6.4 x 


io- 


-40 


7.1 x 


io- 


39 


1.6 x 


io- 


-44 


0.3 


6.1 x 10- 


42 


6.1 x 10 


-42 


2.6 x 


io- 


-40 


5.8 x 


io- 


-40 


7.1 x 


io- 


-40 


6.8 x 


io- 


-40 


8.6 x 


io- 


-40 


7.4 x 


io- 


39 


1.4 x 


io- 


-44 


0.4 


1.6 x 10- 


-41 


7.7 x 10 


-42 


5.6 x 


io- 


-40 


1.3 x 


io- 


-39 


1.1 X 


io- 


-39 


1.0 x 


io- 


-39 


1.1 X 


io- 


-39 


7.6 x 


io- 


39 


1.3 x 


io- 


-44 


0.5 


4.1 x 10- 


41 


9.7 x 10 


-42 


1.2 x 


10 


-39 


3.0 x 


io- 


-39 


1.8 x 


io- 


-39 


1.5 x 


io- 


-39 


1.5 x 


io- 


-39 


7.9 x 


io- 


39 


1.3 x 


io- 


-44 


0.6 


1.1 x 10- 


40 


1.2 x 10 


-41 


2.6 x 


io- 


-39 


6.7 x 


io- 


-39 


2.8 x 


io- 


-39 


2.1 x 


io- 


-39 


2.0 x 


io- 


-39 


8.1 x 


io- 


39 


1.3 x 


io- 


-44 


0.7 


2.7 x 10- 


40 


1.5 x 10 


-41 


5.5 x 


io- 


-39 


1.5 x 


io- 


-38 


4.4 x 


io- 


-39 


3.1 x 


io- 


-39 


2.7 x 


io- 


-39 


8.4 x 


io- 


39 


1.2 x 


io- 


-44 


0.8 


7.1 x 10- 


40 


1.9 x 10 


-41 


1.2 x 


io- 


-38 


3.4 x 


io- 


-38 


6.9 x 


io- 


-39 


4.5 x 


io- 


-39 


3.6 x 


io- 


-39 


8.7 x 


io- 


39 


1.2 x 


io- 


-44 


0.9 


1.8 x 10- 


-39 


2.5 x 10 


-41 


2.5 x 


10 


-38 


7.6 x 


io- 


-38 


1.1 X 


io- 


-38 


6.6 x 


io- 


-39 


4.7 x 


io- 


-39 


8.9 x 


io- 


39 


1.2 x 


io- 


-44 


1.0 


4.8 x 10- 


-39 


3.1 x 10- 


-41 


5.4 x 


io- 


-38 


1.7 x 


io- 


-37 


1.7 x 


io- 


-38 


9.7 x 


io- 


-39 


6.3 x 


io- 


-39 


9.2 x 


io- 


39 


1.2 x 


io- 


-44 



Measured in erg s 1 (M Q yr 1 ) 1 . 
b Measured in erg s~' A~' (M Q yr" 1 )' 
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Table 9. SFS estimators based on A^yc and L1500 for Pop. III. 





IB(2 - 120 Ma) t 


IB(1 - 100 Mo) 


IB(1 - 500 Me) 


Age 


^Lyc" 1 L l500 b 






4 Myr 

5 Myr 

6 Myr 


1.45 x 10~ 47 8.30 xl0~ 34 
1.95 x IO- 47 9.98 xlO- 34 
2.53 x IO- 47 1.14 xlO"" 


1.94 xl0~ 47 1.11 x 10 33 
2.61 x IO- 47 1.34 x 10~" 
3.40 xl0~ 47 1.52 x 10~" 


2.15 xl0~ 47 1.23 xl0- 33 
2.89 x IO- 47 1.48 x 10-" 
3.76 x 10~ 47 1.69 x 10-" 



1 Values for 2-120 M Q were obtained analytically from those for 1 - 100 M ( 

a Measured in s' 1 M G ~'. 

b Measured in erg s _I A -1 M Q -'. 



Table 10. SFR estimators based on Nhyc and L1500 for Pop. III. 





EB(2- 


120 Mo) + 




EB(1 - 


100 M ) 


EB(1 - 


500 Mo) 




Age 


/VLyc" 


j b 

MSOO 




A Ly c a 


j b 

M500 


A'Lyc" 


/ b 

M500 




10 Myr 


1.27 x IO- 54 


7.96 x 10' 


-41 


1.70 x IO- 54 


1.07 x IO- 40 


1.26 x IO- 54 


8.67 x 10" 


-41 


30 Myr 


1.15 x IO- 54 


5.60 x 10" 


-41 


1.55 x IO- 54 


7.50 x IO- 41 


1.18 x IO- 54 


6.61 x 10" 


-41 


250 Myr 


1.14 x IO- 54 


3.53 x 10" 


-41 


1.53 x IO- 54 


4.73 x IO- 41 


1.17 x IO- 54 


4.52 x 10" 


-41 



+ Values for 2-120 M G were obtained analytically from those for 1-100 M Q . 

a Measured in s -1 (M Q yr -1 )-'. 

b Measured in erg s _I A -1 (M Q yr -1 )-'. 

Table 11. IMF correction factors for all SFS estimators, to be applied on values from Tables [T]|4] when IMF mass limits 0.1 - 100 
M or 1 - 100 M are assumed. See Sect.|4]for a more thorough discussion. 



Mass limits 


IB(0.1 - 100 M ) 


IB(1 - 100 M ) 


Corr. factor 


3.41 


1.34 



Table 12. IMF correction factors for SFR estimators, to be applied on values from Tables |5][8] when IMF mass limits of 0.1 - 100 
M or 1-100 M are considered. See Sect.|4]for a more thorough discussion. 



Magnitude 


EB(0.1 - 100 M ) 


EB(1 - 100 Mo) 


^FIR 


3.54 


1.36 


N L yc 


3.94 


1.47 


Ltnd 


3.51 


1.37 


AsoflX 


3.41 


1.34 


-^1500 


3.54 


1.37 


-^2000 


3.53 


1.36 


-^3500 


3.47 


1.34 


-^4400 


3.28 


1.28 


-^5500 


3.26 


1.27 


^-22200 


3.36 


1.32 



The values taken for L softX were calculated analytically. 



Otf-Floranes & Mas-Hesse: Calibration of star formation rate tracers 
Table 13. IMF correction factors for SFS estimators, to be applied on values from Tables [TJ|4] when IMF slope a 







lB(a = 


1) 


IB{a = 3) 


Magnitude 


4 Myr 


5 Myr 


6 Myr 


4 Myr 


5 Myr 


6 Myr 


£-FIR 


0.61 


0.93 


1.22 


2.41 


2.03 


1.81 


N Lyc 


0.54 


0.85 


1.17 


2.76 


2.31 


2.01 


£-Tad 


0.32 


0.52 


0.71 


3.69 


2.99 


2.60 


£-1500 


0.63 


0.97 


1.29 


2.36 


1.99 


1.76 


£2000 


0.61 


0.95 


1.27 


2.40 


1.99 


1.76 


£-3500 


0.55 


0.90 


1.17 


2.49 


1.99 


1.78 


£-4400 


0.50 


0.78 


0.95 


2.55 


2.09 


1.96 


£-5500 


0.49 


0.71 


0.89 


2.59 


2.26 


2.10 


£-22200 


0.48 


0.63 


0.80 


2.75 


2.64 


2.40 



Table 14. IMF correction factors for SFR estimators, to be applied on values from Tables |5][8] when IMF slope a 





EB(a = 1) 


EB(a = 3) 


Magnitude 


10 Myr 


30 Myr 


250 Myr 


10 Myr 


30 Myr 


250 Myr 


£-FIR 


0.37 


0.45 


0.54 


2.55 


2.09 


1.55 


M.yc 


0.25 


0.25 


0.25 


3.68 


3.67 


3.67 


£rad 


0.37 


0.66 


0.74 


2.95 


1.89 


1.75 


£-1500 


0.41 


0.49 


0.60 


2.41 


1.97 


1.49 


£2000 


0.40 


0.50 


0.63 


2.39 


1.92 


1.40 


£•3500 


0.36 


0.46 


0.64 


2.50 


1.92 


1.26 


£<4400 


0.45 


0.64 


1.16 


2.20 


1.61 


0.98 


£5500 


0.46 


0.66 


1.21 


2.27 


1.68 


1.00 


£22200 


0.60 


0.90 


1.52 


2.34 


1.82 


1.16 



